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Abstract 

A simple, seven-parameter trial function is proposed for a description of the ground state of the 
Lithium atom. It includes both spin functions. Inter-electronic distances appear in exponential 
form as well as in a pre-exponential factor, and the necessary energy matrix elements are evaluated 
by numerical integration in the space of the relative coordinates. Encouragingly accurate values of 
the energy and the cusp parameters as well as for some expectation values are obtained. 
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I. INTRODUCTION 



It is well-known that the standard quantum chemistry approaches to calculation of the 
energies of the low-lying states of few-electron atoms are characterized by slow convergence. 
This convergence problem was demonstrated in a discouragingly explicit way for the case of 
the Helium atom in recent studies by Korobov [1] (using a sort of exponential Hylleraas basis) 
and by Schwartz ^ (using various trial functions that included correlated exponentials, 
pre-exponential integer and fractional powers of both nuclear-electron and electron-electron 
distances, as well as logarithmic terms). An unpleasant drawback of these studies is an 
absence of confidence that high accuracy obtained for the energy guarantees a comparable 



accuracy in expectation values 



18| This question is vital when relativistic corrections are 



studied, particularly in view of the fact that some of these corrections are defined by the 
expectation values of singular and/or local quantities (for a discussion see e.g. j^). Recent 
advances in experimental techniques have now led to experimental data whose understanding 
requires an accurate knowledge of the relativistic corrections 5] . As one of possible ways to 
handle this situation, one of the present authors (FEH) has proposed to look for simple, few- 
parameter "ultra-compact" trial functions which guarantee reasonably high overall accuracy. 
One way of characterizing this approach is to describe it as a search for the most accurate 
ew-parameter trial functions. This line of endeavor is illustrated by work on the H2 molecule 
6|, for which a nearly optimum 14-parameter function was reported, and by work on the 



He isoelectronic series Q, where optimum wavefunctions of up to four configurations were 
generated. 

The present contribution deals with a search for an optimum ultra-compact wavefunction 
for the ground state of the Li atom. An important issue for such a study is how to determine 
the overall quality of a trial function. The viewpoint taken in the present work is to use 
as a quality measure the error in the cusp parameters obtained from the trial function 
(residues arising from the Coulomb singularities of the potential). Of course, this criterion 
becomes reasonable only if the cusp conditions are not artificially fixed to their exact values 
by the choice of form for the trial function. We note that the most popular methods for 
atomic computations use Gaussian-type orbitals, and that while such bases can provide 
extremely accurate energies, they usually lead to vanishing cusp parameters and thereby 
have significant drawbacks for the description of relativistic and other local effects. 
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In the work reported here, the trial functions that were examined consist of exponentials 
in all the relative coordinates, in some cases multiplied by pre-exponential factors dependent 
linearly on the interparticle distances. The matrix elements that arise have been evaluated 
numerically by methods used previously by one of the authors 



II. WAVEFUNCTION AND VARIATIONAL METHOD 

The nonrelativistic Hamiltonian for the Lithium atom under the Born-Oppenheimer ap- 
proximation of zero order, i.e. with the Li nucleus assumed to be of infinite mass, is (in 
Hartree atomic units) 

3 „N 3 3 



1=1 ^ ' i=\ j>i ■' 



where Vj is the 3-vector of the momentum of the ith electron, Z is the nuclear charge 
(here Z=3), is the distance between the ith electron and the Li nucleus, and rij are the 
interelectron distances. The kinetic energy part of H, is conveniently written in terms of the 
distance coordinates and Uk = rij,k ^ i ^ j 

y"v2 _y-/^ 2— -— 



i=l i=l 



+ y2(- ^TTl^- ' 2 

^ V UkTi OTiduk 2 UiUk ouiduk J 

as long as the wave function has no explicit angular dependence. The summation P runs 
over 2, j, k which are the six permutations of 1, 2, 3. 

The variational method is used to study the ground state of the Lithium atom. Physical 



relevance arguments are followed to choose the trial function (see, e.g. Turbiner [ll|]). In 
particular, we construct wavefunctions which allow us to reproduce both the Coulomb sin- 
gularities in Ti and in and the correct asymptotic behavior of large distances. As a result 
the wavefunction of the 25*1/2 Li ground state is written in the particular form 

il) = A[(/)(fi,r2,f3)x] , (3) 

with 
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where the pre-exponential factor is a hnear function of its arguments, while and aij are 
(non-hnear) parameters. A is the three-particle antisymmetrizer 

A = / - Pi2 - Pi3 - P23 + P231 + P312 . (5) 

Here Ptj represents the permutation i j and Pijk stands for the permutation of 123 into 
ijk. In Eq. (|3]), x denotes a doublet spin eigenfunction {S = 1/2) expressed as a linear 
combination 

X = X1 + BX2 (6) 

of two linearly independent spin functions spanning the doublet spin space of quantum 
numbers S = 1/2, M, = 1/2: 

Xi = 2-'/^ [a(l)/3(2)a(3) - /3(l)a(2)a(3)], (7) 

and 

X2 = 6"^/^ [2a(l)a(2)/3(3) - /5(l)a(2)a(3) - a(l)/3(2)a(3)]. (8) 

In Eq. iQ, B is a parameter which can be used to obtain the optimum spin function, and 
a{i),f3{i) are spin up, down functions of electron i. In total, the function ip of Eq. is 
characterized by seven parameters, plus any that may occur in the pre-exponential factor /. 

The matrix elements of 7i can be written as integrals over the nine dimensions repre- 
sented by r*!, "'^3- Integrations over three dimensions describing overall orientation are 
easily performed, and we end up with six-dimensional integrals over the relative coordinates 
(ri, r2, r3, ri2, ri3, r23). While it is in principle possible to reduce these integrals analytically 
to one-dimensional integration: to expressions involving dilogarithm functions, as first shown 
by Romm a„d Hill 5. the aualy* properfes of the resultmg expressions were found to 
be quite comphcated (see Harris [13[). For that reason, the primary method used in the 
present research was direct six-dimensional numerical integration. 

These numerical integrations were carried out using a suitable partitioning of the to 
subdomains based on a profile of the integrand (for details, see e.g. Turbiner and Lopez 
In each subdomain the numerical integration is done with a relative accuracy of ~ 10~^ to 
10^^ by use of the adaptive DOIFCF routine from NAG-LIB [l3] in a parallel manner. Due 
to the complicated profiles of the integrands the numerical calculations are very difficult and 
if not done with great care can lead to a serious loss of accuracy. By comparing numerical 
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and analytical evaluations of some of the simpler matrix elements, it was verified that the 
numerical methods were reliable at least to six significant digits. 

Minimization of the energy with respect to the nonlinear parameters was performed using 



the minimization package MINUIT from CERN-LIB 



III. RESULTS 
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The Li ground-state energies obtained for optimized variational wavefunctions of the form 
given in Eqs. Q and (jl]) are displayed in Table [H the corresponding optimized variational 
parameters are given in Table [Tll Each of the first seven rows of Table [T] describes a wave- 
function with a different pre-exponential factor; the last row of the table reports the energy 
obtained from the most accurate existent Li ground-state calculation ^, a result extrapo- 
lated from a 9576-term wavefunction of Hylleraas type and probably accurate to within 10~^ 
a.u. Using these wave functions we calculated the variational energies and also the values 
of the cusp parameters: 



16|, Eq. (18)], which for the exact wavefunction should be equal to —3 when f = fi (the 
electron- nuclear cusp). 

All the wavefunctions described in Table [ITl display the expected electronic shell structure; 
electrons 1 and 2 correspond to a Is^ pair, with an average exponential parameter that 
exhibits only a small degree of screening relative to the bare- nucleus value a = 3. Significant 
energy improvement has been achieved by the "split-shell" description of this electron pair 
(with ai > 3 > a2)- Electron 3 (to zero order the 2s electron) is not optimally described 
by a pure exponential, and great improvement is obtained by giving it a Slater-type orbital 
(STO) description, as in ?/^2; oi even a hydrogenic 2s form, as in ?/^4. Note that the sign of the 
pre-factor parameter /5i produces the node characteristic for a 2s orbital. The data for ip^ 
and ip^ show that inclusion of a linear interelectron distance improves the variational energy 
(as indeed it must), but the improvement is not as striking as that associated with the factor 
r^. Incidentally, most of the improvement associated with the insertion of the factor in 
ipz simply reflects the fact that with high probability, ri3 is similar in magnitude to rs. This 
observation becomes evident when one notes that the ri distribution is Is-like and far more 
localized than is the 2s-like r^, distribution. Further flexibility, as in ipQ (the best three-term 



TABLE I: Li ground-state energy E, cusp parameter CeN [see Eq. ([9])] and the expectation values 
< r^-^ >,< Tij > for the trial function in Eq. ([3]) with various pre- factors /. 





f{ri,r2,r'i,ri2,ri^,r2z) 


E (a.u.) 


— CeN 


< > 


< ^ij > 


^1 


1 


-7.4547 


2.953 


2.1732 


10.0046 






-7.4712 


2.958 


2.1965 


8.9553 


•03 


^13 


-7.4682 


2.955 


2.1922 


8.9457 


^4 


(1 + /9i rs) 


-7.4727 


2.958 


2.2091 


8.6552 




(1+71 ns) 


-7.4686 


2.955 


2.1990 


8.8953 




(1 + /3i rg + 71 ria) 


-7.4451 


2.901 


2.2501 


8.7688 




(1 + ^1 ra + 71 ria) 


-7.4729 


2.961 


2.2027 


8.7330 


'Exact' 




-7.47806^ 


3.0 


2.1982/ 


8.6684-^ 



" /3i = -2.44486 

71 = 9.53316 
" Pi = -3.82716, 71 = -0.47333 

Pi = -2.7T713, 71 -0.26645 
^ from Ref. 3| (rounded) 

from Ref. (rounded). 

prefactor) gives little additional gain over that already achieved in 1IJ4. All the wavefunctions 
also exhibit small negative values of the interelectron screening parameters aij, thereby 
improving the description of the repulsive electron-electron correlation. However, neglecting 
this screening by setting all aij = (see the function ipl, Table II) worsens significantly the 
variational energy as well as the cusp parameter, see Table I. Note that a function of this 
type was used in [5], yielding the most accurate variational energy thus far obtained, but at 
the expense of a very long expansion. 

All the results reported in this paper use the doublet spin function that optimizes the 
trial energy for the given spatial function. From the small values of B in Table II, we see 
that in every case the dominant contributor to the spin state is (as expected) that which 
couples the Is and Is' spatial functions into a spin singlet. However, inclusion of the other 
member of the spin basis does influence the energy to some extent; for ipi (the trial function 
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TABLE II: Variational parameters and Uij in [a.u.] ^ and B (dimensionless) for some trial 
functions from Table I. 







1p2 


■04 


re 






3.2892 


3.3065 


3.3038 


3.3051 


3.3044 


a2 


2.3343 


2.3291 


2.3519 


2.0657 


2.3603 


as 


0.4336 


0.7004 


0.7473 


0.6690 


0.7327 


ai2 


-0.2108 


-0.2050 


-0.2150 


0. 


-0.2218 


ai3 


-0.0411 


-0.0311 


-0.0194 


0. 


-0.0227 




-0.0404 


-0.0316 


-0.0313 


0. 


-0.0276 


B 


0.06295 


0.01416 


-0.00201 


0.00277 


-0.00202 



with the largest optimum B), use of the second spin state lowers the trial energy by about 
0.001 a.u. 

Recently, it was shown by one of the authors that a correct treatment of the domain of 
WKB asymptotics of the wavefunction at large distances is very important for getting a high 
quality trial function. Usually, the large-distance asymptotic expansion of the exponential 
phase of the wavefunction contains several terms that grow as the distance increases. All 
these terms should be reproduced in a trial function, otherwise exponential deviation from 
the exact function at large distances occurs. In the case of Lithium it can be verified that 

if = -log^! = aifs + a2logr3 + 0(1), ^ oo , r 1,2 fixed , (10) 

where ai^2 are constants. These two terms in the expansion are the only terms that grow 
as r3 increases, and failure to reproduce them in the trial function can lead to an expo- 
nentially large deviation of the trial function from the exact eigenfunction at large (and 
intermediate) distances r^. All six of the functions we study reproduce at least the linear 
term in Eq. (fTOl) and all but ipi reproduce both terms. This observation might be considered 
as an explanation why ipi gives worse results for energy than the other wavefunctions (see 
Table [T]) . Unfortunately, at this time a complete analysis of the asymptotic behavior of the 
exponential phase at large distances (in different directions in 6D r-space) is missing; such 
an analysis would be helpful for the identification of adequate trial functions. 
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Following the above analysis, it is not surprising that the function ipi gives the largest 
deviation from the exact values in both the energy (~ 0.5%) and the cusp parameter CeN 
(~ 1.6%), while for the most accurate function ipQ these deviations are respectively ~ 0.01% 
and ~ 1.3%. It is worth noting that for the functions ipi-Q an increase of the accuracy in 
energy corresponds to a decrease of the error in the cusp parameter CeN (see Table [I]). A 
similar situation occurs for the expectation value (vij) in comparison with the value reported 
in jj]. For {r~j^) the largest deviation from the result from \^ occurs for the function ipi 
which also provides the largest deviation for energy, cusp parameter and < >. We 
must emphasize that all deviations in energy as well as in the expectation values occur 
systematically at the third significant digit. 



IV. CONCLUSION 



Simple and compact few-parameter trial functions are presented for the ground state 
of the Lithium atom. These already provide a very accurate ground state energy. These 
functions are the most accurate among existent few-parameter trial functions. However, 
the presented analysis does not seem final since one could explore the more extended pre- 
f act or 

/(ri, ra, rg, r^, ri3, rsa) = {I + pi ri + P2 ^2 + /^a '"a + 7i ^23 + 72 + 73 ris) , 

which seems beyond the computer resources presently available to us. 

The wavefunctions used in the present work can be easily modified to study excited states 
of the Lithium atom. 
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